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1. INTRODUCTION 

In recent years, there has been important research activity in the fault-tolerant control for Hybrid 
systems. These systems contain continuous and discrete systems described by differential equations, 
automata and sequential logics. To control these systems even in case of failure, we adopted the fault 
tolerance control (FTC). The FTC is developed in two approaches: Passive approaches based on a robust 
controller that ensures the imperceptibility to the systems, and active approaches that change control 
operations for each detected fault [1], [2], [3]. 

In this paper, we are interested in a new control strategy for piecewise affine systems (PWA) that is 
used for system's reconfiguration. The contribution in this paper aims essentially to reconfigure the non linear 
hybrid system using the Hammerstein method that works only for linear systems. The linearization is 
obtained by using Taylor method. 

The next sections in this paper are organized as follows: Section 2 presents the stability concepts, 
Section 3 describes the linearization, Section 4 discusses the reconfiguration with the Hammerstein method 
and Section 5 presents the simulation example. 


2. STABILITY NOTIONS 

In this part, we will present the stability technique with lyapunov method. 
Lyapunov method 

To verify the stability of an hybrid system using the lyapunov method or more generally the LMI 
method. [3][4]. There exists a matrix 


P=P'>0 C1) 
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We say that the system is quadratically incrementally stable if it verifies the following LMI equation. [5] [6]. 


PA +A'P<0 ,i=1,..., p (2) 


3. LINEARIZATION (TAYLOR METHOD) 
If a Non Linear system has multiple dynamic variables on discrete time, it will be written 
as the following: [7], [8], [9], [10], [11], [12]. 


x fi(Xi -Xash ot) 
= f(x,u)=] : (3) 


x fa (Xi -Xps Uje- Um) 


where (X,...X,,) is the state vector and (U,...U,,,) is the control vector; n and m are positive integers. 


When the system operates around an equilibrium point and the signals involved are small, then it is 
possible to approximate the Non Linear system with Taylor series by computing the Jacobian matrices in 
order to obtain a linear system with the following state space model: 


X= Ax + Bu 
| (4) 


y = Cx + Du 


Where Ais the state matrix, Bis the control input matrix, C is the output matrix, Dis the feedthrough 


matrix and y is the output vector. Now, we will approximate f (x, u) with Taylor series: 


j(u-u,) (5) 


Xie 
Xe 5|: is a vector determined by solving: 
Xhe 
fi(Xi -Xas ee ES) Gaid e 
, (6) 
fa (Xie Xps Ue Un) ui = ne 
To check the stability, it is necessary to find the Jacobian matrix: 
Of(%u) Of, (xu) 
; OX, OX, 
X,U ‘ 
Ok u) e 
Ox 


fu) Af Cu) 
OX, Ox 


n 
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that verifies: 


det(J, —AI) =0 (8) 
Of (x,u 
Where: J, = of (xu) ; 
OX (Xeste) 
A, 
=|: is the eigenvector, and [ is the identity matrix. 
A 


n 


The equilibrium is locally stable if (—1<A<1),. 
Finally, we determine y = g(x,u) 


0g (x,u 0g(x,u (9) 
y=g(x,,u,)+ e lea) è Jl atzu) 
pe ay) eee 
ay [lve ou (tree) 
(10) 
C _ 29(%u) D _ 6g(xu) 
Ox |(Kette)” Gu |(Xette) 


4. RECONFIGURATION METHOD 
In this part, the piecewise affine system which is a class of a hybrid system is defined. Moreover, 
the reconfiguration method applied to hybrid system is processed. 


4.1. Piecewise Affine System (PWA) 

The piecewise affine system is characterized by dividing the state-space in a finite number of 
regions and associating to each one an affine linear equation. [13]. The nominal system is modeled by the 
following PWA model: 


ee = Ax(t)+b, + Bu,(t)+ B; d(t) 
(11) 
y(t) = Cx) 
Where: d (t): disturbance 
Ba : disturbance Input matrix 
bi : The affine terms 
A; .the real matrices of appropriate dimensions for all i. 
In this section the reconfiguration method is processed. 
4.2. Fault Model 
The fault event suddenly changes the nominal PWA system to the faulty PWA system: 
x, (t)= Ax, (t)+B, u,(t)+ B; d()+b,, > 


yŒ) = C,x;(t) 
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Where: B, = corresponds to faulty actuator 


f 
C r= corresponds to faulty sensor 


U; = corresponds to faulty controller 


In case of actuator faults, we change the measurements of the matrix B; as follows: 

B, =(1- Beta) x B (13) 
In case of Sensor faults, we change the measurements of matrix C pas follows: 

C, =(alpha)x C (14) 


Beta, alpha are random values that keep the system faulty in the same interval of operation as the nominal 
system. 


4.3. Reconfiguration with Hammerstein Method 

The Hammerstein method is modeled by new input and output matrices of a system having actuator 
faults or sensor faults [14]. The reconfiguration with Hammerstein method is defined by the following 
equation: 
1. In case of actuator faults, we obtain the following virtual actuator PWA: 


3O = A (0) +b, + By, u, (t) 
y(t) =C X(t) (15) 
u; (t) =M x(t) + Nx(t) + B; b, 


With 

B,, =+(B-B,M) (16) 
X is found by the following equation: 

(A -B,M)+(A-B,M)'x <0 (17) 


2. in case of a sensor fault, we obtain the following virtual sensor PWA: 


R (t) = A,R, (t) +b, + Bu,(t)+ Ly, (t) 


(18) 
Asi = (4 —LC;) 
With X is found by the following equation: 
2, (A -C,L)+(A-C,L)'%, <0 (19) 


We find M and L with the following steps. 
There exist matrices X,Y and X,Y, that gratify the LMIs: 
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X, =X] >20 
VEY 20 
(20) 
XA+tAX, -Y C; -C Y7 <0 
(i =1,..., p) 
and 
X, =X] 20 
Ter 
T T T < (21) 
AX,+X,A —B,Y,-Y, B; <0 
(i =1,..., p) 
M is a perturbation observer for the defective system: (for actuator fault) 
LO X `Y, (22) 


The Hammerstein method is summarized in the following algorithm. 
Steps 1 to 3 describe the nominal system before failure, in step 4 the defects are detected, the gains is 
calculated in step 5, in step 6 to 8 the reconfigured closed-loop system is run. 


Hammerstein Method 
1) request : PWA model Ą,, b,, B, C 


2) begin the nominal closed-loop system 


C; =C, B; =B,b;;=b,,L=0,M =0 
Ni = 1,x(0)= x 


3) execution the nominal closed-loop system up to actuator or 
sensor fault f detected and isolated 


4) make the fault model b, ,B,,C, and update the PWA (12) 


5) find LMI (20),(21) and calculate (22) M O Y, X „~^ for 


actuator fault, (23) LU X a Y, for sensor fault 


6) calculate By ; (16), X (17) for actuator fault, X (19) for 


sensor fault 


7) Update PWA actuator fault (15), PWA sensor fault (18) 


8) execution reconfigured closed-loop system 


5. EXPLICATIVE EXAMPLE 
An application of Algorithm Hammerstein to the model of a two tanks with 4 modes is presented in 
this section. 


5.1. The Non Linear Hybrid System 

In this part, we will present a model of two tanks with 4 modes as shown in Figure 1. The two-tank 
system is a typical example of a non linear hybrid system. It has been considered as a reference problem for 
the diagnosis and detection of failures in [15]. Table 1 shows the parameters of system: 
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Table 1. List of Symbols 
Symbol Declaration 

Tı Tank 1. 

To Tank 2. 

hy liquid level in tank 1, initially 0.5m 

hy liquid level in tank 2, initially empty 

P Pump controlled in all or none. 

Gh Cp, C3, Cy Conducts. 

Vi Valve kept open during the functioning system 

V2 Valve kept open during the functioning system. 

Va Controlled Valve (V4is opened at t = 240s and closed at t = 380 s). 

qp Debit in the pump P, (qp= 0 when homax= 0.2m and qp = 0.001m°/s when homin= 0.1m). 

qı Debit in valve V4. 

q2 Debit in valve V>. 

q3 Transfer debit of conduct C3 at 0.5 m height. 

qa Debit in valve V4. 

S Section area of Tı and Tz; S=0.0154 m°. 

A Conducts Section; A= 0.000036 m’. 


Figure 1. Two tanks system 


The functioning of this system is described in Figure 2 by the automate hybrid. 


. 1 . 1 
h=<5@-%-%) k= 5 Qp- 4-2-4) 


. 1 . 1 
h, =—(q, h, =—(q,+4:) 
h 5h) 2 z 13 


h)<0.5, hy<0.5, V4 Closed h,>0.5, ho<0.5, V4 Closed 


1 


ghetii) 


h)<0.5, ho<0.5. V4 Opened h:>0.5, ho<0.5, V4 Opened 


ex Opening V4 
ex Closing V4 


Figure 2. Automate hybrid 
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The debits q4, q2, q3, q, are given by the following equation: 


q,=A/2gh, 
q2 =sign(h,-h, )A,/ 2g |h,-h,| 
q,=sign(h,-0.5)A./2g|h,-0.5| 


q,=A./2gh, 


(24) 


Note that the system is discretized with a sampling period T, = 0.5s. 
Figure 3 represents the liquid level variation of the non linear System. 


0 200 400 600 800 1000 


Figure 3. Liquid level variation of the non linear system 


In this paper, we are interested in the reconfiguration of the non linear hybrid system. Whereas, the 
Hammerstein method is valid only for linear systems. We can summarized the work as: (1) providing an 
linearization method for a non linear Hybrid Systems ;(2) developing a Hammerstein algorithm which 
addresses both of the actuator faults and the sensor faults. 


5.2. Linearization of the non Linear Hybrid System 
This method of linearization is applied to our hybrid system: 


ah, Oh oh, oh, 
a| l | % os 
oh, oh, ôq, Od, 


1 0 0 0 h 
C } o-[ } a a u=] 
0 1 0 0 h, q, 
Matrices searching process for each mode: 


1. Mode 1: h, <0.5, h, < 0.5 , V, Closed 
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The system's equations will be written as follows: 


- 1 
h, = 54 — 4-H) 


i (26) 
h, = z 
- 1 
h, = (q, — AŅ2gh, —sgn(h, —h, )AJ2g |h, z h, |) 
see i 
h, = senth, —h,)A 2g |h, — h,|) 
The steady state operating points h, and h,, are given by solving: 
h, (hehe) = he 
(27) 
h, (hehe) = he 
The Jacobian matrix: 
ah êh, E OEE OEE | 
Fz oh, oh, E A29 Jh, „lh, > h, | lh, z h,| 
dh, dh 2s 


1 1 
dh, Oh, vlh, -4,| lh, —h,| 
The Jacobian matrix at the operating points h,, and h,,: 
oh, oh, 
dh, |(Metee) ap, | (herhre) es 0.1776 
J = 55 


| a, ah, p 
ch, (heshe) ch, 


0.1776 -0.1776 


(he „hoe) 


Calculating the eigenvector of J, : 


yE G l E a. 
A j (0.4471 


(—1< 4 <1) and (—1< A, <1) so the equilibrium is locally stable, then A = J, . 
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The control matrix B is calculated as follows: 


oh, oh, 
1 
oq, 4, = 0 
B=| . . |=|s 
ch, ah, 6 
ôq p ôq p 


2. Mode 2: h, > 0.5, h, < 0.5 , V, Closed 


The system's equations will be written as follows: 


>. 1 
ae tp q-q, -q;) 


. 1 
h, ee + qs) 


> 1 
h,=—(q,-A [2gh, ~sgn(h, —h,)A,/2g |h, — h,| -sgn(h, —0.5)A,/2g |h, —0.5|) 
- 1 

h, =—(sgn(h, —h,)A, 2g |h, -h,|+sgn(h, -0.5)A,/2g|h, -0.5]) 


The Jacobian matrix: 


1 1 1 


J _Ay2g Jh hh, -h,| ih, -0.5| 
j 2s 1 n 1 


The Jacobian matrix at the operating points: 


—0.1259 
J 
0.0555 


a: 
la, = h, | 
1 


j lh, -h,| 


0.0481 
—0.0481 


Calculating the eigenvector of J, : 
—0.0223 
a,) (0.1517 
(—1< 4 <1), so the equilibrium is locally stable, then A = J, . 


The control matrix B is calculated as follows: 


w 
II 
oO újļ=e 
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(28) 
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3. Mode 3: h, >0.5, h, <0.5 , V, Opened 


The system's equations will be written as follows: 


= oj 
h= q-q, -q;) 


. 1 (29) 
h, = sh +q; =q) 


- 1 
h, = Vi — A,/2gh, —sgn(h, —h,)A,/2g |h, — h,| —sgn(h, —0.5)A,/2g |h, — 0.5]) 
- 1 

h, =< (senth, —h,)A,/2g |h, — h,| + sgn(h, —0.5)A,/2g |h, —0.5| —q,) 


The Jacobian matrix: 


1 1 1 1 
2s 1 


a Oe est +s 
la, —h,| Jh, P 0.5| hh, = h,| 


The Jacobian matrix at the operating points h,, and h,,: 


e 


—0.1191 0.0399 
= 
0.0473 -—0.0399 


Calculating the eigenvector of J,: 
—0.0207 
A p= A; ce 
A, ) \ 0.1383 
(—1< A <1), so the equilibrium is locally stable, then A= J,. 


The control matrix B is calculated as follows: 


° 
| 
| 


4. Mode 4: h, < 0.5 , h, < 0.5 , V, Opened 


The system's equations will be written as follows: 


- 1 
h = 5 —U-) 
(30) 


. 1 
h, = 5h =) 
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uoi 
h, =— (q, — AV2gh, -sgn(h, — h,) AJ 2g |h, — h|) 
eA 

h, = Cen —h,)AŅ2g |h, —h,| — da) 


The Jacobian matrix: 
ee ee, aa oO 
ayaa} Jh vin-rl vik -nl 
2s 1 z 1 
lh, —h,| ah, -h| 


The Jacobian matrix at the operating points: 


ee a 
Jos 


J= 


0.0681 —-0.0681 


Calculating the eigenvector of J,: 


0.0446 
~ | -0.2654 


(—1< 4 <1), so the equilibrium is locally stable, then A= J,. 


Ne 
I 
ON 
Gees 


The control matrix B is calculated as follows: 


° 
| 
| 


5.3. Reconfiguration method 

The adopted reconfiguration method applied to hybrid system is processed in this section. In the first 
part, an actuator fault is introduced to the system and the performance of the Hammerstein method by 
reconfiguration of the actuator fault is checked. In the second part, a sensor fault is inserted to the system and 
then a reconfiguration of a sensor fault is undertaken. 


5.3.1.Actuator Fault and Reconfiguration: 


a. Actuator Fault model: 
The new matrix Br equals: 


B, =(1- Beta)x B 
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Figure 4 represents the linear system and the actuator fault system. 


Linearized System VS Actuator Fault System 


weeeeeee H1 Linear System 
— =H2: Linear System 
== H1: Linear Actuator Fault 
~~" H2: Linear Actuator Fault 


y 


Figure 4. System with actuator fault 


Figure 5 is Figure 4 focused on the fault part of the actuator fault system. 


Linearized System VS Actuator Fault System 
H HE pan sy pn H1: Linear System 
= =H2: Linear System 
H1: Linear Actuator Fault 
—— He: Linear Actuator Fault 


06 


05 


o4 


+ + + 4 + + + + + 


450 500 550 600 650 700 750 800 850 900 950 1000 


Figure 5. System with actuator fault (Zoom) 


Failures detection is neatly visible on the curves of Figures 4 and 5, the fault is introduced during the time 
interval [400-1000s], the failure is detected at t=460s, t=660s, and t=982s. 

b. Reconfiguration of the Actuator Fault by Hammerstein method: 

We find M (perturbation observer): 


eg [0.4214 -0.5342 
MOY,X,7=10 


0.1500 -0.5717 


X Is found by the following equation: 


%(A-B,M)+(A-B,M)'X <0 
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~ 0.0069 0.1053 
x= 
0.1053 0.1676 


Figure 6 represents the reconfiguration for an actuator fault system. 


Actuator Fault System VS Reconfigured System 
~~" H1: Actuator Fault System 
H2: Actuator Fault System 
—— H1: Reconfigured System 
~ H2: Reconfigured System 


Figure 6. Reconfiguration for an actuator fault 


Figure 7 shows the Comparison of the linear system with the reconfiguration system for an actuator fault. 


Linear System VS reconfigured Actuator Fault System 
H — =H1: Linear System 
2 s H2: Linear System 
——— H1: Reconfigured System 


Figure 7. Comparison of the linear system with the reconfiguration system for an actuator fault 


Figure 8 is a zoomed version of the Figure 7; it shows that the adopted method has successfully compensate 


the effect of the faults on the system. 
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Linear System VS reconfigured Actuator Fault System 


¿J= =H1: Linear System 
| oss H2 Linear System 
——— H1: Reconfigured System 
~ H2: Reconfigured System 


400 500 600 700 800 900 1000 


Figure 8. Comparison of the linear system with the reconfiguration system for an actuator fault (Zoom) 


Figure 9 shows the efficiency of the reconfiguration with Hammerstein method for an actuator fault. 


Actuator Fault Error 
0.0114435 
0.011443 
0.0114425 
0.011442 
0.0114415 
0 500 1000 0 500 1000 


Figure 9. Error between the linear system and the reconfiguration system for an actuator fault 


The computation of the error is realized for a time period between the linear system and the 
reconfiguration system for an actuator fault. It is interesting to note that the error’s value is close to zero. 
From the previous results, the figures show that the Hammerstein method is very effective to reconfigured 
the actuator faults. 


5.3.2. Sensor Fault and Reconfiguration: 


a. Sensor Fault model: 
The new matrix C; equals: 


1 0 
C, =(alpha)<C= alpha 5 i] 
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Figure 10 represents the linear system and the sensor fault system 


Linear System VS sensor Fault t System 
F t + e H1: Linear System 
] |= =H2: Linear System 
— H1: fault System 


Figure 10. System with sensor fault 


Figure 11 is a zoomed version of the Figure 10. 


= ©H1: Linear System 
= =H2: Linear System 
— H1. fault System 
— H2: fault System 


ee a a a 


400 420 440 460 480 500 520 540 560 580 600 


Figure 11. System with sensor fault (Zoom) 
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Failures detection is neatly visible on the figures 10 and 11, the fault is introduced during the time interval 
[400-1000s], the failure is detected at t=440s, t=460s, t=475s, t=490s, and t=533s. 
b. Reconfiguration of the Sensor Fault by Hammerstein method: 

We find L (state observer) 


L 


X Y, =-0.2274 
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With X is found by the following equation: 


R (A -C,L)+(A-C,L)'2, <0 


sa 102 0.8099 0.0108 
X= 
0.0108 0.8126 


Figure 12 represents the reconfiguration for a sensor fault system. 


fault sensor System VS Reconfigured System 


— H1: fault sensor System 
— =H2: full sensor System 
——— H1: Reconfigured System 
—— H2: Reconfigured System 


Figure 12. Reconfiguration for a sensor fault 
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Figure 13 represents the comparison between the linear system and the reconfiguration system for a sensor 


fault. 


linear System VS Reconfigured System 

===- H1: linear System 

sere" H2: linear System 

— H1: Reconfigured System 


Figure 13. Comparison of the linear system with the reconfiguration system for a sensor fault 
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Figure 14 is a Zoomed version of the igure 13; it shows that the effect of the faults on the system is 
compensated by the use of the adopted method. 


linear System VS Reconfigured System 


cocco H1: linear System 
se... H2: linear System 
——— H1: Reconfigured System 
~ H2: Reconfigured System | 


0.15} 


Figure 14. Comparison of the linear system with the reconfiguration system for a sensor fault (Zoom) 


faults. 


Figure 15 represents the error between the linear system and the reconfiguration system for a sensor 
fault; it shows the efficiency of the reconfiguration with Hammerstein method for the sensor fault system. 


Sensor Fault Error 


10 


(ae) 0% ae 


Figure 15. Error between the linear system and the reconfiguration system for a sensor fault 


The computation of the error is realized for a time period between the linear system and the 
reconfiguration system for a sensor fault. It is interesting to note that the value of the error is around zero. 
From the previous figures, the Hammerstein method gives a very effective reconfiguration of the sensor 


Comparison between the reconfiguration methods for non linear hybrid systems and other methods: 
to prove the efficiency of the proposed reconfiguration method, an evaluation is given in this section. 
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The result obtained using the Hammerstein method is compared to the methods developed in [16], [17] and 
summarized in the Table 2. 


Table 2. Comparison of Methods 


Method Hammerstein Method in [16], [17] 

Stability For each method, we try to compute the stability using from Lyapunov method 
(Matlab LMI Control Toolbox). The methods guarantee the stabilily of the system. 

Time for detected failure Fast (T= few seconds) Slow 

Reconfiguration the effect of the faults on the system is minimizes the propagation of failure 
compensated effects 


The previous results show the effectiveness of the Hammerstein method for non linear hybrid 
systems. This approach ensuring a very important result to reconfiguration method (for the actuator faults and 
the sensor faults) and ensuring the stability of the systems. 


6. CONCLUSION 

During last years, there has been a lot researchs in the active fault tolerant control for non linear 
systems and linear Hybrid systems. [18], [19], [20], [21]. The contribution in this work is essentially to 
reformulate a new approach to improve the reconfiguration performance of the non linear hybrid system. 
This approch compensate the effect of faults and ensure the stability of the closed-loop system. 

Therefore, this paper developed the linearization of a non linear Hybrid Systems using the Taylor 
method. Moreover, the active fault tolerant control based on the Hammerstein method has been proven to be 
efficient for the reconfiguration of the non linear Hybrid Systems in case of the intervention of the actuator 
fault and the sensor fault. The application of this method on the hydraulic system (two tanks with 4 modes) 
gives hopeful results.The future research will be focused on the reconfiguration method using the fault 
tolerant control for non linear hybrid systems based on neural network. 
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